Current-voltage characteristics of quasi-one-dimensional superconductors: An S-curve 

in the constant voltage regime 
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Applying a constant voltage to superconducting nanowires we find that its IV-characteristic ex- 
hibits an unusual S-behavior. This behavior is the direct consequence of the dynamics of the 
superconducting condensate and of the existence of two different critical currents: j C 2 at which the 
pure superconducting state becomes unstable and j c \ < j C 2 at which the phase slip state is realized 
in the system. 
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The majority of the experiments on the resistive state 
in quasi-one dimensional systems were performed in the 
constant current regime and at temperatures close to T c . 
It is extremely difficult to apply voltage to a supercon- 
ductor because the current density induced by the ap- 
plied electric field inevitably reaches the critical value 
and destroys superconductivity in the sample. The de- 
crease of the superconductingcurrent by the appearance 
of phase slip centers Q, 0, EJ [J is not effective in this 
case because of the large heating of the sample at low 
temperatures. At temperatures close to T c the heating 
can be suppressed due to the low value of the critical cur- 
rents but in this case the applied electric field does not 
penetrate deep into the sample because of the existence 
of regions near the N-S boundary where the drop of the 
applied voltage occurs 0, . 

This situation drastically changes with the appearance 
of nano-technology and the ability to create long (to al- 
low the appearance of phase slip centers) superconduct- 
ing wires with a small cross section (to decrease the effect 
of heating) . In this Letter we present results on the be- 
havior of such nanowires in the constant voltage regime. 
We found that the I-V characterestic in this case has 
a remarkable S-shape. Our theoretical analysis based on 
the time-dependent Ginzburg-Landau equations (TDGL) 
shows that such a behavior is a direct consequence of 
the dynamics of the superconducting condensate and we 
predict new unusual features which still need additional 
experimental study. 

The superconducting nanowires were prepared by elec- 
trodeposition into nanopores of homemade track-etched 
polycarbonate membranes 0- For the lead nanowires, a 
22 /im thick membrane (with pore diameter ~ 40 nm 
and pore density ~ 4-10 cm -2 ) and an aqueous solution 
of 40.4 g/1 Pb(BF 4 ) 2 , 33.6 g/1 HBF 4 and 15 g/1 H 3 B0 3 
were used 0, while in the case of the tin nanowires, a 50 
fj,m thick membrane (with pore diameter ~ 55 nm and 
pore density ~ 2T0 9 cm -2 ) and an electrolyte of 41.8 g/1 
Sn(BF4)2 in water solution were applied. Constant po- 



tential of -0.5 V versus an Ag/AgCl reference electrode 
was used in a three-electrode configuration in order to re- 
duce the Pb 2+ or Sn 2+ ions into the nanopores. As shown 
in Fig. 1 , the nanowires are cylindrical and the diameter 
is uniform along their length. In order to perform elec- 




FIG. 1: SEM-FEG micrograph of lead nanowires after 
dissolution of the hosting polycarbonate membrane with 
dichloromethane. 

trical transport measurements on a single nanowire (still 
kept inside the membrane), a self-contacting technique 
was developed. In addition of the thick gold cathode (~ 
1 /im) from which the nanowires start to grow up inside 
the membrane, another thin gold layer (in the range of 
50 to 200 nm) was deposited on the other side exposed 
to the electrolyte prior to the electrodeposition. This 
contacting layer covers the surface but in contrast to the 
cathode, it is not thick enough to plug the pores. As 
the nanowires do not grow at the same speed, the first 
emerging nanowire interrupts the growth of the others by 
favoring the formation of a film of the deposited material 
on this contacting layer. 

All the transport measurements are done in an electro- 
magnetically shielded environment and special precau- 
tions were taken to reduce the noise. The voltage across 
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FIG. 2: Current-voltage characteristics at different tempera- 
tures of a) a Pb nanowire (diameter 40 nm, length 22 /im) and 
b) a Sn nanowire (diameter 55 nm, length 50 (jm). Results 
are shown for the current driven mode (black curves) and the 
voltage driven mode (grey curves). 



the sample was measured by a Keithley voltmeter. From 
the numerous measurements performed both on Pb and 
on Sn nanowires, it appears that almost no depression of 
T c was observed. In addition, a superconducting transi- 
tion was found to occur when the residual resistivity ratio 
(RRR) is larger than about 2 while for smaller RRR val- 
ues an insulating behavior is observed 0| . Our nanowires 
are moderately disordered and the RRR of the samples 
was larger than 10. Structural characterization of the 
Pb nanowires using TEM experiments confirm also their 
polycrystalline character [lfj . 

In Fig. 2, the DC current-voltage (I-V) characteris- 
tics of the Pb and Sn nanowires are reported, applying 
either a DC current or a DC voltage. All the curves pre- 
sented in Fig. 2 were measured in both directions and no 
hysteresis was observed. These measurements were also 
performed with reversed polarity without any changes 
compared to the reported results. In the usual current 
driven experiment we observe for the Pb nanowire (Fig. 
2(a)) the successive appearance of two PSCs at low tem- 
perature. Applying the phenomenological SBT model 0, 
the size of the first PSC can be estimated to be about 
18 /J,m, which is twice the quasiparticles diffusion length. 
In spite that successive PCSs tend to avoid those already 
in place, the second PSC created in this Pb nanowire is 
thus forced to interpenetrate the first one, which explains 



why the second jump in resistance is smaller. The cur- 
rent driven experiment on the 50 fj,m long Sn nanowire 
(Fig. 2(b)) also shows the formation of a PSC at low 
temperature but for this particular material, the PSC 
extension is larger (around 40 fim) and the formation 
of a second PSC almost coincides with the transition to 
the normal state. Interesting new features were observed 
when measuring the reverse way, i.e. applying constant 
voltage and measuring the current. Here, the current 
flowing into the nanowire was determined by measuring 
the voltage across a 1 f2 resistance added in series and 
the voltage developed across the sample was measured 
separately. In this voltage driven experiment, we found 
an interesting S-shape behavior which occurs at low tem- 
perature in the formation region of these PSCs, both for 
the Pb and the Sn nanowires. 

To understand the experimental results we used the 
generalized TDGL equation Q which describes the dy- 
namics of the superconducting condensate if) = \ip\e 10 
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(i) 

The physical quantities are presented in dimcnsionlcss 
units (see e.g. Ref. Q) and we neglected the vector 
potential, because in our case the self-induced magnetic 
field is negligible in size. Eq. (1) is supplemented with 
an equation for the electrostatic potential tp which is ob- 
tained from the condition for the conservation of the to- 
tal current in the wire, i.e. divj = 0. The parameter 7 
depends on the material and we took U ~ 5.79 Q. 

Although Eq. (1) is valid only close to T c (see Q) 
we expect to find clues of the dynamics of the system 
even far from T c . This expectation is motivated by the 
fact that results found from the solution of the stationary 
Ginzburg-Landau equations for mesoscopic superconduc- 
tors [TlT | agree with experiment even far from T c . 

If we apply an external current to the wire then the 
system exhibits hysteretic behavior 0, Hj . If one starts 
from the superconducting state and increases the current 
the superconducting state switches to the resistive super- 
conducting or normal state at the upper critical current 
jc2- With decreasing current it is possible to keep the 
sample in this state even for currents up to j c \ < j C 2 
(which we call the lower critical current). The supercon- 
ducting resistive state is realized as a periodic oscillation 
of the order parameter in time at one point of the su- 
perconductor 0, ■ When the order parameter reaches 
zero in this point, a phase slip of 27r occurs. Such a state 
is now called a phase slip state and the corresponding 
point a phase slip center (PSC). The meaning of the cur- 
rent jd is that for j < j c \ the phase slip solution cannot 
be realized (in the absence of fluctuations) and thus the 
current j c \ is the critical current at which PSCs start to 
appear. Thus there is a limiting cycle (see Ref. 01) f° r 
Eq. (1) giving rise to this process. 
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FIG. 3: The current voltage characteristics of a supercon- 
ducting wire of length L — 40£ and 7 = 10 (Aq ~ 2.3£). 
Dotted(solid) curve for the I(V) = const regime. 

The existence of such two critical currents is, in the 
V = const regime, responsible for a complicated I-V 
characteristic (see Fig. 3). The I(V) curve has an overall 
S-like behavior which at low voltages exhibits an oscilla- 
tory behavior superimposed on it. We found that these 
unusual properties are typical for superconducting wires 
with 7 > 1. To understand the physics leading to this 
unusual behavior we consider the dynamics of ip in the 
wire for different voltages. 

When we apply a voltage V to the wire of length L a 
constant electrical field appears in the sample E = jg = 
V/L. As a result the superconducting condensate will 
be accelerated and the momentum p — V# will increase 
in time. For small voltages and electric fields the ab- 
solute value of the order parameter and current density 
are approximately described by quasi-equilibrium expres- 
sions with \tp\(t) = l-p(t) 2 and j s (t) = p(t)(l - p(t) 2 ) 
(p(t) ~ j E t) while p < p c (at I/( > 1 we have 
p c ~ 1/V3). For p > p c the above spatial independent 
solution becomes 'unstable' (like for the case of a ring 
in the presence of an applied magnetic field as discussed 
111 Refs. EH). The order parameter vanishes in the 
center of the wire and a jump of the phase by 2tt occurs 
in this point. It decreases p by Ap = 2tt/L (and hence 
the superconducting current density) in the entire wire. 
A second phase slip will appear in the same point and 
the phase slip process repeats itself as long as j > j c \. 
When the current density becomes less than j c \ the or- 
der parameter will increase till its quasi-equilibrium value 
\4>\ (t) = 1 — p(t) 2 with p(t) < p c . Than the applied elec- 
tric field will again accelerate the superconducting con- 
densate and the above process is repeated (see Fig. 4). 
We define the time which is necessary to decrease the or- 
der parameter from \ip\ — -\/2/3 till zero and again back 
to the value |^>| = \j2fi as the transition time T tr and 
call this entire process the transition process. 



FIG. 4: Dependence of the order parameter in the center of 
the wire (in the minimal point of |i/j(s)|) as function of time. 
Solid(dotted) curve corresponds to the voltage in point 1(2)- 
V = 0.0024 (V = 0.0028) of Fig. 3. 



It is easy to show that at small voltages (in the limit 
V — > 0) the number of phase slips which occurs dur- 
ing the transition process in an ideal wire is given by the 
following expression N m j n = Nint ((p c — p c i)(L/27r + 1)), 
where p c \ is the smallest real root of equation j c i — 
p c i(l—Pci) and Nint(x) returns the nearest integer value. 
The number of phase slips occurring during the transition 
process increases with increasing wire length. 

Because of the presence of an external electric field, 
p also increases during the transition process. Roughly 
speaking, if the voltage in the wire is so large that 
jETtr > 2tt/L the number of phase slips during the tran- 
sition will increase by one (see Fig. 4) . It leads to an in- 
stant decrease of the average current. The reason is that 
the maximum current during current oscillations jmax 
for a small change of V does not change but the mini- 
mum value of the current jmin decreases by about 2tt/L, 
and hence < I > should decrease. This effect is most 
pronounced in short wires where N m in is small and the 
term 2w/L is large. If N does not change with increas- 
ing V then < I > increases because both minimum and 
maximum currents grow. 

When the period T = 2nN/V becomes of order 2T tr 
the dynamics of the system changes considerably (at volt- 
age V\ in Fig. 3). Starting from this point the average 
current < I > monotonically decreases with increasing 
applied voltage. For these values of the voltages the sys- 
tem is unable to return to the "quasi-equilibrium" state 
(it is a state such that if we switch off the voltage in- 
stantly the system remains in this state for an infinitely 
long time in the absence of fluctuations). At any mo- 
ment of time the order parameter is more suppressed in 
the point where the oscillations occur. At these volt- 
ages jmax decreases with increasing V and in principle 
it should reach the current j c \ (and j m in should reach 
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FIG. 5: Time dependence of the order parameter in the cur- 
rent decreasing regime for a long wire (L = 40£). Dot- 
ted(solid) curve corresponds to the voltage in point 4(5)- 
V = 0.08 (V = 0.2) of Fig. 3. 



j c i from the bottom). But this is only possible in an 
infinitely long wire. In a finite length wire every phase 
slip event leads to a decrease of the momentum p and 
the current j by a finite value. As a result the lower 
critical current in the V = const regime is equal to 
jYi — jci+dj (with dj ~ 1/i) at which the 'usual', like in 
the / = const regime, periodic phase slip processes with 
period T = 2tt/ < V > starts. This corresponds to the 
voltage Vi in Fig. 3. 

In long wires a complex behavior is found in the current 
decreasing region due to the existence of several phase 
slip events during the transition process. At the initial 
parts of this region an irregular behavior is typical. In 
Fig. 5 we show (dotted curve) an example of such a 
situation (at V = 0.08 - point 4 in Fig. 3). Not a single, 
but several periods T = 2t:N/V are present. This regime 
is replaced at higher voltages (for example at V = 0.2 - 
point 5 in Fig. 3) by a regime with a single period and 
the situation becomes similar to the case at low voltages. 



We expect that the S shape of the TV characteristic in 
the V = const regime will be conserved even for temper- 
atures far from T c , because the only condition for its oc- 
currence is the existence of two different critical currents 
jci < jc2- While V < Vi the current in the wire will be 
an oscillating function with maximal amplitude j C 2 and 
minimal amplitude j c \. It means that the time averaged 
current should be larger than j c \. At V > Vi we expect 
that jmax decreases due to the same mechanisms as was 
discussed above and hence the average current should de- 
crease. But the maximal current cannot decrease below 
jd and for V > V% both j m m and j ma x have to increase 
with increasing voltage. 

In our theoretical consideration we neglected the effect 
of thermo- and/or quantum-activated phase-slip centers. 



We believe that this is the reason why in the experi- 
ment no hysteresis was observed in the / = const regime. 
The larger these fluctuations the larger they will decrease 
the measured j%£ p and increase the measur ed (see 
also the discussion of this question in Ref. and in 

principle they may suppress the hysteresis. But in the 
V = const regime the situation is quite different. If the 
time needed for the appearance of a fluctuating phase 
slip center is much larger than the period of oscillation 
of the current then we can neglect their effect. 

In conclusion, we presented the first experiment on the 
I-V characteristic in the voltage constant regime for nar- 
row superconducting Pb and Sn wires. An unexpected 
S-shape response is found which is absent in a current 
driven measurements. A theoretical analysis based on 
the generalized TDGL equations gives a qualitative ex- 
planation of this phenomena which is driven by a highly 
nonlinear time response of the superconductor to elec- 
tric fields leading to the periodic creation of PS centers. 
Furthermore, we predict an oscillatory I-V characteristic 
at low voltages which was not found in the experiment 
presumedly due to the large length of the wire and hence 
the relatively large value of N m i n . 
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